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The internal thermal noise in LIGO's test masses is analyzed by a new technique, a direct appli- 
cation of the Fluctuation-Dissipation Theorem to LIGO's readout observable, x(t) = (longitudinal 
position of test- mass face, weighted by laser beam's Gaussian profile). Previous analyses, which 
relied on a normal-mode decomposition of the test-mass motion, were valid only if the dissipation 
is uniformally distributed over the test-mass interior, and they converged reliably to a final answer 
only when the beam size was a non-negligible fraction of the test-mass cross section. This paper's 
direct analysis, by contrast, can handle inhomogeneous dissipation and arbitrary beam sizes. In 
the domain of validity of the previous analysis, the two methods give the same answer for S x (f), 
the spectral density of thermal noise, to within expected accuracy. The new analysis predicts that 
thermal noise due to dissipation concentrated in the test mass's front face (e.g. due to mirror coat- 
ing) scales as by contrast with homogeneous dissipation, which scales as 1/ro (ro is the beam 

on : radius); so surface dissipation could become significant for small beam sizes. 
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I. INTRODUCTION 



' Random thermal fluctuations are expected to be the 
dominant noise source for the first interferometers in 
the Laser Interferometer Gravitational Wave Observa- 
tory (LIGO) at frequencies between 35 and 100 Hz. This 
thermal noise is generally decomposed into a suspension 
thermal noise and an internal thermal noise for the test 
masses. The former can be traced back to the friction in 
the test masses' pendular suspension system; the latter 
is due to internal damping inside the test masses thcm- 
ON ■ selves. Traditionally, themal noise calculations have been 
\ based on a normal-mode expansion Q], [pj . However, 
Q^- Gonzalez and Saulson have also performed an exact cal- 
culation of the suspension thermal noise by applying di- 
6JT). rectly the Fluctuation-Dissipation (FD) theorem || in 
it's most general form, due to H. B. Callan and T. A. 
Welton U . The purpose of this paper is to use the gen- 
eral method of Gonzales and Saulson to calculate the 
J-j | internal thermal noise. 

In Section II we analyze a general situation when a 
measuring device (e.g. a laser interferometer) monitors 
the displacement of the surface of a test mass whose in- 
ternal degrees of freedom are in themal equilibrium with 
each other. We develop a general formalism for using the 
FD theorem to calculate the thermal noise in the most 
general surface readout quantity. In brief our method is 
as follows: 

To work out the thermal noise at a particular frequency 
/, one should mentally apply pressure oscillating at this 
frequency to the observed surface of the test mass. The 
spatial variation of this pressure should mimic that of 
the light beam intensity (for example, in the case of a 
gaussian beam this oscillating pressure has a gaussian 
profile of the same widthi as the beam). The thermal 
noise is then given by 

2k B T Wmss m 
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where ks and T are the Boltzmann's constant and the 
temperature of the mirror respectively, Fq is the ampli- 
tude of the oscillating force applied to the surface (i.e. 
the pressure integrated over the surface), and Wdiss is 
the time-averaged power dissipated in the test mass when 
this oscillating pressure is applied. 

To demonstrate the computational power of this gen- 
eral approach, in Section III we consider the case of a 
cylindrical fused silica test mass monitored by a circular 
gaussian laser beam. For the case when the radius of 
the beam is much less then the size of the test mass and 
the dissipation is uniformly distributed throughout test 
mass volume, we derive an analytical expression for the 
thermal noise [cf . Eq. ([l5]) of Section III ] : 



Here c, Eq, and <j> are the Poisson ratio, Young's mod- 
ulus, and dissipational loss angle [Eq. (|ll|)1 of the test- 
mass material, r$ is the radius of the laser beam (which is 
defined here as a radius at which the intensity of light is 
1 1 e of the maximum intensity) , R is a characteristic size 
of the test mass, and / = 1.87322... in the case of a gaus- 
sian beam. Putting numbers in Eqs. (|l|) and (^|), we find 
that our results are in agreement with those of Raab and 
Gillespie j3|, who used the more complicated and com- 
putationally involved method of normal-mode decompo- 
sition. It is interesting to note that as tq/R tends to 
zero, our simple analytical formula becomes more pre- 
cise, whereas the more complicated and computationally 
involved method of normal-mode decomposition requires 
summing over a larger number of modes and thus be- 
comes computationally more expensive. 

Not only can the normal-mode decomposition be com- 
putationally expensive, it can also be misleading. We 
demonstrate this point in Section IV by considering a 
test mass which has a lossy surface, e.g. due to a lossy 
mirror coating. We estimate the contribution of the sur- 
face to the thermal noise using the general method of 
Section II, and show that it differs from the estimate ob- 
tained by the method of normal modes (which gives a 
result too small by a factor of at least ~ tq/R). This 
breakdown of the normal-mode analysis will in general 
happen when the sources of friction are not distributed 
homogeneously over the test mass. The fundamental rea- 
son is that in this case different normal modes can have 
a common Langevin driving force (which is not so if the 
defects are distributed homogeneously). 

Our analysis shows that thermal noise due to surface 
losses near the laser beam spot scales as S x (f) oc \jr\, 
whereas thermal noise due to volume losses scales as 1/tq. 
Correspondingly, for small beam spots the surface losses 
could become significant. To protect against this, it is 
important to keep the surface near the laser beam spot 
as free of potential sources of friction as possible. 
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II. GENERAL METHOD 



For concreteness, consider a situation where LIGO's 
laser beam is shining on the circular surface of one of 
LIGO's cylindrical test masses. The phase shift of the 
reflected light contains information about the motion of 
the test mass's surface. The variable read out by this 
procedure can be written as 

x(t) = J f(r)y(f,t)d 2 r. (3) 

Here r is the transverse location of a point on the 
test-mass surface, and y(r, t) is the displacement of the 
boundary along the direction of the laser beam at point 
f and time t. The form factor f(r) depends on the laser 
beam profile and is proportional to the laser light inten- 
sity at the point r 0; it is normalized by J f(r)d 2 r = 1. 

The internal thermal noise of the test mass is defined 
as the fluctuations in x(t), and our objective is to find the 
spectral density S x (f) of these fluctuations. Wc assume 
that the test mass is in thermal equilibrium at tempera- 
ture T. 

Callen and Welton's general- 

ized Fluctuation-Dissipation Theorem || says that the 
spectral density of the fluctuations of LIGO's readout 
variable x(t) is given by the formula 

S x {f)= 1 %\Re[Z{f)}\, (4) 

where ks is Boltzman's constant and Z(f) is a com- 
plex impedance associated with x(t). This complex 
impedance can be understood and computed as follows. 
Introduce a special set of generalized coordinates for the 
test mass's degrees of freedom-a set for which x is one 
of the coordinates. (Since x is not the coordinate of a 
normal mode of the test mass, these generalized coordi- 
nates will not be the usual ones associated with normal 
modes.) Apply to the test mass a generalized force F(t) 
that drives the generalized momentum conjugate to x 
but does not drive any of the other generalized momenta. 
This generalized force will show up as the following in- 
teraction term in the test mass's Hamiltonian: 

Hint = -F(t)x. (5) 

This driving force, together with the test mass's inter- 
nal elastic forces and internal dissipation, will generate a 
time evolution x(t) of the observable x. Denote by F(f) 
and x(f) the Fourier transforms of the (arbitrary) driv- 
ing force F(t) and the observable's response x(t). Then 
the impedance that appears in the thermal noise formula 
Eq. d) is 

Z(f) = 2mfx(f)/F(f). (6) 

The physical nature of the driving force F(t) can be 
deduced by inserting the definition (0) of the observable 
x into the interaction Hamiltonian (pf): 
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where 



Hmt = - J P(7)y(r,t)d 2 r, (7) 



P{f,t) = F(t)m. (8) 



From Eq. (0) we see that the generalized force F(t) con- 
sists of a pressure P(f,t) [Eq. (§h] applied to the test 
mass's surface. Note that the spatial distribution of this 
pressure is the same as LIGO's laser beam intensity pro- 
file. 

The real part of the impedance, Re[Z(f)], describes 
the coupling of the test mass's dissipation to the observ- 
able x. We can see this most clearly by applying an os- 
cillatory pressure P{r,t) = F cos(2ir ft) f(f) to the test 
mass's face. From the response formula ([|) we infer that 
the power Wdiss that this oscillatory pressure feeds into 
the test mass, and that the test mass then dissipates, is 
related to \Re[Z(f)]\ by 

|Jfc[Z(/)]| = ^=. (9) 

Substituting Eq. (g) into Eq. (g), we get 



2k B T W di! 
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Equation (|l^) is the most important equation of this 
paper. Let us reephasize it's physical content: 

1. Apply an oscillatory pressure 

P(r,t) = F cos(2n ft)f(r) to the face of the test 
mass. 

2. Work out the average power Wdiss dissipated in the 
test mass under the action of this oscillatory pres- 
sure. 



3. Use Fq and Wdiss in Eq. (10) to calculate S x (f) 



This procedure is different from the one employed in 
previous calculations of internal thermal noise for the 
LIGO and VIRGO test masses @, §, g|. The previ- 
ous authors decomposed a test mass's motion into nor- 
mal elastic modes; then they calculated the contribution 
of each mode to S x independently and added up these 
contributions. This method of "normal-mode decompo- 
sition" works fine in many cases, but it has two draw- 
backs: 



1. The fundamental assumption in this method is that 
different normal modes have independent Langevin 
forces. This assumption is correct only if the 
sources of friction are homogeneously distributed 
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over the test-mass volume. It breaks down if the 
defects are more concentrated in one place than 
in others — for example, when there is significant 
damping concentrated in the test-mass surface. We 
will return to this in Section IV. 

2. For a small laser beam diameter the sum over nor- 
mal modes converges very slowly, so one has to sum 
over many modes, which may be computationally 
expensive. By contrast, using the new method de- 
scribed in this paper, one can write down a simple 
analytic expression for the low-frequency noise in 
the case of a narrow laser beam. In the next sec- 
tion we derive this expression and make compar- 
ison with the normal-mode decomposition results 
derived in S. 

III. THERMAL NOISE DUE TO 
HOMOGENEOUSLY DISTRIBUTED DAMPING 

Consider the case where all the friction in the test mass 
comes from homogeneously distributed damping. It is 
conventional to characterize such friction by an imagi- 
nary part of the material's Young's modulus: 

E = E [1 + i<f>(f)] ; (11) 

</>(/) is called the material's "loss angle". It is sus- 
pected H , ||] that for fused silica, which will be used in 
LIGO's test masses, <fr might be independent of frequency 
within LIGO's detection band (but there is no evidence 
for such behavior of cf> for figh-quality resonators — see 
for some healthy scepticism) . In this /-independent case 
the damping is called "structural" . 

To calculate the thermal noise for homogeneous dissi- 
pation, we express Wdiss in Eq. ( |lC)| ) as 

Wd iss = 2tt/IW<K/), (12) 

where t/ max is the energy of elastic deformation at a mo- 
ment when the test mass is maximally contracted or ex- 
tended under the action of the oscillatory pressure of Eq. 

®- 

LIGO's detection frequencies (10 — 300Hz) are much 
lower than the eigenfrequencies of the test mass's normal 
modes (the lowest of which is ~ 6kHz) ; so we can assume 
constant, non-oscillating pressure P(r) — Fof(r) when 
evaluating t/ max . 

In the case when the beam profile is gaussian and the 
centre of the light spot coincides with the centre of the 
transverse coordinates, we have 

f(r) = -^e-^^o, (13) 

where tq is the radius of the laser beam. When the char- 
acteristic size of the test mass R is much greater than 
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ro, we can approximate the test mass as an infinite half- 
space in order to find t/ max . Appendix A uses elasticity 
theory to derive C/ max in this case [cf. Eq. flAq)]: 



F, 



U m ^ = ^—{l-a 2 )I 1 + 0(^1 , (14) 



n 2 E r 



R 



where Eq and a are the Young's modulus and Poisson 
ratio of the material respectively, and I ~ 1.87322. Here 
O(ro /R) is a correction due to the finite size of the cylin- 
der. Putting Eqs. @ and (|l|) into Eq. ©, one gets 



Below we take the numerical values [] used by Gillespie 
andRaab §|: r = 1.56cm, E = 7.18xl0 10 Pa, a = 0.16, 
<f> = 10~ 7 , the mirror diameter of 25cm and the mirror 
length of 10cm. Gillespie and Raab, after summing over 
the relevant ~ 30 modes, get 

5^ R (100Hz) ~ 8.0 x 10- 40 m 2 /Hz. (16) 



Our analytical approximation ([L5|) (which should be valid 
to within ~ 10 percent in this case) gives 

S^lOOHz) ~ 8.7 x 10" 40 m 2 /Hz. (17) 

Notice that our analytic expression in Eq. ( |l5| ) gets more 
exact when tq/R — > 0, whereas, by contrast, the sum over 
modes converges more slowly and gets more complicated. 

The ratio tq/R may turn out to be of order unity in 
real experiments. In this case, Eq. ( |l5] ) can only be used 
for order-of-magnitude estimates. To work out the exact 
value of the internal thermal noise, one would need to 
calculate U max numerically. We have done such a numer- 
ical computetion using finite-element techniques. More 
specifically, we have used finite-element software called 
PDEasc2D [Version 3.0], which runs as part of Mascyma 
[Version 2.1], to solve the elasticity equations for the 
loaded mirror and to compute C/ max and, by virtue of 
Eqs. ( |l2| ) and (JT^) , S x . The exact answer for the mirror 
and light spot parameters given above is 

S^lOOHz) = 8.76 x 10" 40 m 2 /Hz, (18) 

which is consistent (better than expected) with our ana- 
lytical approximation. 

The purpose of the present section is to convince the 
reader that the method presented in this paper is correct 
and could be computationally cheaper than the normal- 
mode expansion. The next section concentrates on the 
cases where a direct application of the FD theorem can 
be crucial for getting the right results, and the method 
of normal-mode decomposition fails. 



1 Note that our definition of the beam radius (location where 
intensity has fallen to 1/e of its central value) differs by y/2 
from the beam radius of Ref. Q (location of 1/e amplitude 
falloff). 
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IV. THE CASE OF SURFACE DAMPING 

In this section we study thermal noise due to surface 
losses-caused, e.g., by inadequate polishing or by a lossy 
mirror coating. 

From Eq. ( [10| ) we see that the key quantity in the ther- 
mal noise calculation is the power dissipated in the test 
mass when an oscillating pressure is applied to the laser 
beam spot on the test-mass surface. The power dissi- 
pated at each point of the material is proportional to 
the square of the stress at this point. Most of the sur- 
face stress is in or near the spot to which the pressure is 
applied, so 

V 'o / 'o 

Thus the thermal noise due to the surface damping scales 
like 

S x (boundary) tx 1/rg. (20) 

For comparison, the thermal noise due to bulk damping 
[Eq. (|l5|)] scales as 

^(bulk) oc l/r . (21) 

Thus as the spot size decreases, the thermal noise due 
to surface damping grows faster than that due to bulk 
damping. 

Contrast this conclusion with the intuition one gets 
from normal-mode decomposition. There one is con- 
cerned with how much the surface contributes to the 
quality factors (Q's) of the normal modes. For a typ- 
ical mode the strain at the surface is at most of the 
same order as the characteristic strain inside the test 
mass (likely, much less for first few modes — because of 
the free boundary condition). Therefore, one would pre- 
sume that the surface contributes no more than some 
mode-independent fraction of the test mass's Q's. In or- 
der of magnitude this fraction should be the ratio of the 
power dissipated in the surface to that in the bulk if one 
applies an oscillating pressure uniformally to the whole 
surface, which in the context of our method corresponds 
to a beam radius of R. Therefore the normal-mode es- 
timate of the surface thermal noise is at least ro/R less 
than the correct value. 

Current experiments show that the mirror coating does 
not contribute significantly to the Q's of the test-mass 
normal modes. The conclusion commonly made is that 
coating is also not likely to contribute significantly to the 
internal thermal noise. The above analysis shows that 
this conclusion is not justified and that there might be 
a significant contribution of the coating to the internal 
thermal noise, despite the fact that Q's are not signifi- 
cantly changed. 
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V. DISCUSSION AND CONCLUSION 

The normal-mode decomposition of the thermal noise 
is exact when the defects are distributed homogeneously 
through the volume of the test mass. However, as was 
shown explicitly in Section IV for the case of surface 
losses, when the defect distribution is not homogeneous, 
the normal-mode decomposition may be misleading, and 
a direct application of the Fluctuation-Dissipation theo- 
rem is required. 

Thermal noise is ultimately linked to friction in the test 
mass; this friction is caused by various (structural and 
otherwise) defects. Those defects which are closer to the 
beam spot will contribute more to the thermal noise that 
is read out by the laser beam's phase shift. Although this 
fact is a direct consequence of the formalism developed 
in this paper, we would like to give an intuitive example 
in order to emphasize this point. 

Consider, for the sake of simplicity, a one dimensional 
elastic test mass with two identical defects A and B, 
as shown on Fig. 1; A is closer to the beam spot than 
B. Each of these defects creates a random stress which 
pushes apart or pulls together the left and right (relative 
to the defect) parts of the test mass. By conservation of 
momentum, the part of the test mass which is lighter will 
respond more to the random stress than the other part; 
therefore defect A will have a larger effect on the optical 
readout than the B. 

Note that if the defects A and B are positioned sym- 
metrically with respect to the centre of the test mass, 
they will have the same effect on the Q's of all elas- 
tic modes (we assume for simplicity that only one- 
dimensional longitudinal modes are present — and all of 
them are either symmetric or antisymmetric with respect 
to the centre). Therefore, the normal-mode decomposi- 
tion applied to the test-mass with just one defect — A 
or B — would give the same result for the thermal noise 
as read by the laser. Clearly, we have found yet another 
illustration of the breakdown of the normal-mode decom- 
position . 

The considerations presented above lead to the follow- 
ing advice for real experiments: keep the neighbourhood 
of the laser beam spot as clean of defects as possible. 

Not only does our direct applicatin of the Fluctuation- 
Dissipation Theorem have broader validity than the 
normal-mode decomposition; it is also be computation- 
ally simpler. In the case of homogeneous structural 
damping it yields a simple analytical expression for the 
internal thermal noise spectrum [cf. Eq. (J15|)] : 



This result is consistent with the numerical sum-over- 
modes done in Ref. J3| and is accurate when the radius 
of the laser beam is small relative to the size of the test 
mass, i.e. in the regime when the sum over modes con- 
verges especially slowly. When ro/R is not small, a nu- 




(22) 
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merical solution of the elasticity equations to deduce the 
dissipation power Wdiss! & n d thence the thermal noise 
(10), is straightforward and is probably also much sim- 
pler than performing a sum over modes. 
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APPENDIX A: THE STRAIN ENERGY IN A 
TEST MASS SUBJECTED TO A GAUSSIANLY 
DISTRIBUTED SURFACE PRESSURE 

The objective of this Appendix is to derive Eq. (|lj) of 
Section III for the energy of elastic strain in a cylindrical 
test mass when the pressure P(r) = Fof(r) is applied to 
one of it's circular faces. (As was discussed in Section 
III, we can assume that the pressure is constant in time 
since LIGO's detection frequencies are much lower than 
the lowest normal- mode frequency). For a circular laser 
beam with a gaussian intensity profile f(r) is given by 
[cf. Eq. ©] 



m = -Ke-^i, (Al) 

'" n 



where we assume that the centre of the light spot coin- 
cides with the centre of the test mass's circular face. 

If the radius of the laser beam r$ is small compared 
to the size of the test mass, we can approximate the test 
mass by an infinite clastic half-space. Then our calcu- 
lation of the elastic energy is correct up to a fractional 
accuracy of O(ro/R), where R is the characteristic size 
of the test mass. 

Let y(r) be the normal displacement of the surface at 
location r under the action of the pressure P(f). In the 
linear approximation of small strains 

y(r) = [ G{r 7 f>)P(fi)d 2 r', (A2) 



where G(r,r') is a Green's function. The calculation of 
G is a non-trivial albeit standard exercise in elasticity 
theory ]^], which gives 

G(f,P)= 1 ~f 1 (A3) 
7r£/o \ r — r \ 

where a is the Poisson ratio and Eq the Young's modulus 
of the material. The elastic energy stored in the material 
is 
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C/max = / P(f)y{r)d 2 r 

1 - a 2 f P(r)P{?) 
ttE J \f-p\ 

I -a 2 2 



TT EqTq 



d 2 rd 2 r' 



-(r +r )/r 



(A4) 



-,d 2 rd 2 r' 



2rr' cos 



where is the angle between f and r' . The integral in 
the last term of Eq. (A.4) (as was pointed out by Glenn 



Sobermann) can be taken by introducing "polar" coor- 



dinates R and 



r = Rcosi 



= R sin d>. One then 



integrates out the radial part of the integrand and ex- 
pands the remaining angular part in a power series with 
respect to cos#; termwise integration of this power se- 
ries finally yields Eq. p4j ) [up to a fractional error of 
0{r o /R)] 



ir 2 E a r 



(l-a 2 )I, 



where 



7T 3 / 2 



E 



(4n- 1)!! 



' (2n)!4"(2n + 1) 



(A5) 



1.87322. (A6) 



It can be shown that if, instead of an infinite half-space, 
we consider a finite cylindrical test mass, the leading 
fractional correction to the elastic energy is of the order 
O(r /R). 
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FIG. 1. Identical defects A and B create fluctuating strees 
in different parts of the test mass. The stress created by defect 
A will influence the phase shift of the laser beam readout more 
than the stress created by defect B, although both A and B 
make identical contributions to Q's of the test mass's elastic 
modes. 
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